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Abstract 

This paper studies near-critical evolution of the spherically symmetric scalar 
field configurations close to the continuously self-similar solution. Using an- 
alytic perturbative methods, it is shown that a generic growing perturbation 
departs from the Roberts solution in a universal way. We argue that in the 
course of its evolution, initial continuous self-similarity of the background is 
broken into discrete self-similarity with echoing period A = ^/2■K = AAA, 
reproducing the symmetries of the critical Choptuik solution. 
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I. INTRODUCTION 



Critical phenomena in gravitational collapse have been a relatively recent and interesting 
development in the established field of general relativity. Following the numerical work of 
Choptuik on the spherically symmetric collapse of the minimally coupled massless scalar 
field [0, critical behavior was discovered in most common matter models encountered in 
general relativity, including pure gravity null fiuid [§] and, more generally, perfect fiuid 
as well as more exotic models. 

The essence of critical phenomena in general relativity is the fact that just at the thresh- 
old of black hole formation, the dynamics of the field evolution becomes relatively simple 
and, in some important aspects, universal, despite the complicated and highly non-linear 
form of the equations of motion. In analogy with second order transitions in condensed 
matter physics, the mass of the black hole produced in near-critical gravitational collapse 
scales as a power lawQ 

MbM oc \p-p*\^, (1) 

with parameter p describing initial data, and the mass-scaling exponent (3 is dependent only 
on the matter model, but not on the initial data family. The critical solution, separating 
solutions with black hole formed in the collapse from the ones without a black hole, also 
depends on the matter model only, and serves as an intermediate attractor in the phase space 
of solutions. It often has an additional symmetry called self-similarity, in either continuous 
or discrete fiavors. 

Discovery of critical phenomena in gravitational collapse was the first real success of 
numerical relativity in which a physical effect was observed in simulations without being 
first predicted by theoretical physicists. For the theoretician, the challenge and attraction 
of studying critical phenomena lies in the possibility of exploring a new class of exact solu- 
tions of Einstein's equations, having simple properties and high symmetry, but previously 
undiscussed. Another interesting thing about critical solutions is that they are very relevant 
to the cosmic censorship conjecture, the long-unsolved problem of general relativity. With 
their ability to produce arbitrarily small black holes and, in the critical limit, curvature 
singularity without an event horizon, in the course of quite generic gravitational collapse, 
they may serve as an acceptable counterexample to the cosmic censorship conjecture (see 
101 and references therein). 

Universality of the near-critical behavior has been explained by perturbation analysis and 
renormalization group ideas IQ-Q, and is rooted in the fact that critical solutions generally 
have only one unstable perturbation mode. In the course of evolution of the near-critical 
initial field configuration, all the perturbations modes contained in it decay, forgetting details 
of the initial data and bringing the solution closer to critical, except the single growing mode 
which will eventually drive the solution to black hole formation or dispersal. In this sense, the 
critical solution acts as an intermediate attractor in the phase space of all field configurations. 



Usually, but there are models with mass gap in black hole production, most notably Yang-Mills 
field, whose behavior is more analogous to a first order phase transition. 
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Because there is only single growing mode, the codimension of the attractor is one. The 
eigenvalue of the growing mode determines how rapidly the solutions will eventually depart 
from critical, and it can be used to calculate the mass-scaling exponent (3. 

As we have mentioned, the critical solution often has additional symmetry besides the 
usual spherical symmetry, called continuous or discrete self-similarity. This symmetry es- 
sentially amounts to the solution being independent of (in case of continuous self-similarity) 
or periodic in (in case of discrete self-similarity) one of the coordinates, a scale. The role of 
this symmetry in critical collapse is not understood at all. Some attempts at finding critical 
solutions made a continuously self-similar ansatz and hit a jackpot 0], whereas others stud- 
ied discretely self-similar solutions from phenomenological point of view . Yet the simple 
question of why a particular matter model should have this or that version of self-similarity 
incorporated in the critical solution still remains a mystery. 

This paper attempts to shed some light on the subject by investigating the dynamics 
of formation of discretely self-similar structure in the gravitational collapse of a minimally 
coupled massless scalar field. As a base point of our investigation, we consider a certain 
continuously self-similar solution, known as the Roberts solution, as a toy-model of the 
critical solution in the gravitational collapse of a scalar field. This solution was constructed 
as a counterexample to the cosmic censorship conjecture [^] and was later rediscovered in 
the context of critical gravitational collapse [0,^. While not a proper attractor [12 



this simple solution resembles in some of its properties more complicated critical solutions 
known only numerically. The aim of the present work is to show how, at least in the linear 
approximation, the discretely self-similar structure arises dynamically in the scalar field 
collapse. The advantage of our approach is that, due to the simple form of the Roberts 
solution, calculations can be carried out analytically and so provide additional independent 
insights different from numerical treatments. 

Using linear perturbation analysis and Green's function techniques, we study evolution 
of the spherically symmetric scalar field configurations close to the continuously self-similar 
solution. Approximating late-time evolution via the method of stationary phase, we find 
that a generic growing perturbation departs from the Roberts solution in a universal way. 
In the course of the evolution, initial continuous self-similarity of the background is broken 
into discrete self-similarity by the growing perturbation mode, reproducing the symmetries 
of the Choptuik solution. We are able to calculate the echoing period of the formed discretely 
self-similar structure analytically, and its value is close to the result of numerical simulations. 



II. THE ROBERTS SOLUTION 

The starting point of our investigation is the Roberts solution, which will serve as a 
background for linear perturbation analysis. It is a solution describing gravitational collapse 
of a minimally coupled massless scalar field, described by the Einstein-scalar field equations 

R,u = 2<P,,<j),., (2) 
□0 = 0, (3) 

which is spherically symmetric and also continuously self-similar. The latter symmetry 
means that there exists a vector field ^ such that 
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££,gpiy = '^QiJ.u, £^4> = 0, (4) 

where £ denotes the Lie derivative. Under these assumptions the field equations can be 
solved analytically, which is most easily done in null coordinates [|l^,|ll|,0. Self-similar 
solutions form a one-parameter family. As the parameter is varied, spacetimes both with 
and without a black hole occur. The solution just at the threshold of black hole formation^] 
is given by the metric 



-2 du dv + dQ^ 



where 



V u'^ — UV, (j) 



In 



1 - 



u 



(5) 



(6) 



The global structure of the corresponding spacetime is shown in Fig. |1|. The infiux of the 
scalar field is turned on at the advanced time f = 0, so that the Roberts spacetime is 
smoothly matched to Minkowskian fiat spacetime to the past of this surface. The junction 
conditions, required for continuity of the solution there, serve as boundary conditions for 
the field equations (|]) and (H). More detailed discussion of this important point is provided 
in Appendix 0. 

The evolution of perturbations of the Roberts solution is most easily followed in a coordi- 
nate system exploiting scale-invariance of the background, so that the self-similarity becomes 
apparent. Therefore, we introduce new coordinates, which we will call scaling coordinates, 
by 

^lnK--l, s = -ln(-u), (7) 



X 



with the inverse transformation 



u 



u 



-^fe^" - 11 



(8) 

The signs are chosen to make the arguments of the logarithm positive in the region of interest 
(w > 0, M < 0), where the field evolution occurs. In these coordinates the metric (^) becomes 



g^^ dx^dx"" = 2e2(^-') [(1 - e~'^'')ds'^ - 2dsdx] + dQ 



and the Roberts solution (Rl) is simply 



X. 



(9) 



(10) 



Observe that the scalar field (j) does not depend on the scale variable s at all, and the only 
dependence of the metric coefficients on the scale is through the conformal factor e~^*. This 
is a direct expression of the geometric requirement in scaling coordinates; the homothetic 
Killing vector ^ is simply — 



^In the early works [ p^ , pA| term critical has been used to designate this solution. This is somewhat 
confusing because, strictly speaking, this solution is not an intermediate attractor of codimension 
one, and so is not critical in the usual sense. Perhaps the term threshold gives better description 
of its nature. In any case, since we are not concerned with the other solutions from the self-similar 
family in this paper, we will refer to the solution (|5|,^) by the name "the Roberts solution" . 
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III. GAUGE-INVARIANT PERTURBATIONS OF THE ROBERTS SOLUTION 



Since we are ultimately concerned with the dynamics of the breaking of the fields away 
from the Roberts solution, the effect due to the growing perturbation modes, we will only 
consider spherically symmetric perturbations here. Non-spherically symmetric perturbations 
decay and so do not play a role in the critical behavior. In this section, we outline how 
spherically symmetric perturbations of the Roberts solution are described in gauge- 
invariant formalism. A general spherically-symmetric metric perturbation is 

Sg^ydx^dx^ = kuu du^ + 2kuv du dv + k^^ dv^ + r^K dVl? ^ (11) 

while general perturbation of the scalar field is 

5(t) = if. (12) 

Under a (spherically-symmetric) gauge transformation generated by the vector 

e = {A,B,0,0), (13) 

the metric and scalar field perturbations transform as 

Ag^i, = £^gf,u, = £^(j). (14) 

The explicit expressions for the change in the perturbation amplitudes under the gauge 
transformation generated by the vector ^ are 



Ah — 


-2B. 




Aky^Tj = 


~^,« 




At — 


-2A, 




r^Air = 


(2m- 


-v)A- 


2r2A(^ = 


vA - 


uB. 



(15) 

Out of four metric and one matter perturbation amplitudes one can build the total of three 
gauge-invariant quantities, one describing matter perturbations 

K If 

/ = y - + ^ y k^^dv, (16) 
and the other two describing metric perturbations 

p = [r -\- kyy kuv uk^u^^jl -\- (2it v^kyy^u/2, (1"^) 

= kuv ~ 2 y ^^^'^ '^^ ~ 2 I 0-^) 
The linearized Einstein-scalar field equations 

=40(,^50,,), 5(n0) = O (19) 
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can then be rewritten completely in terms of these gauge-invariant quantities. It is possible 
to show that the field equations reduce to one master differential equation for the scalar 
field perturbation, 

2u{u - v)Uv + {2u - v)f,, - uf,u - 2/ = 0, (20) 

and two trivial equations relating metric perturbations to the scalar field perturbation 

= 2/,„ + 2f/u, p = 0. (21) 

Once the gauge-invariant quantities are identified, one is free to switch between various 
gauges. We conclude this section by discussing two particularly convenient choices. 

Field gauge {K = k^^ = 0): The scalar field perturbation coincides with the gauge- 
invariant quantity / in this gauge, and expressions for other gauge-invariant quantities 
simplify considerably: 

/ = -<^, 

P ku^ uky^iiu / 2y 

c" — kuv 2 ^ ^uu,v du. (22) 

The linearized Einstein-scalar field equations are at their simplest in this gauge, and the 
derivation of the master equation for / above is almost transparent. The metric and scalar 
field perturbation amplitudes are trivial to obtain: 

^ = -/, 

kuv 2y, 

kuu,v = -4:f/u. (23) 

Null gauge (fcutt — kyy — 0): This gauge was used in the original analysis of spherically- 
symmetric perturbations of the Roberts solution [0. The motivation behind this gauge 



choice is that coordinates u and v remain null in the perturbed spacetime. The expressions 
for gauge-invariant quantities are quite simple here as well: 

f = K/2- ^, 

P (^ -^),™ ky^^, 

a = kuv. (24) 

For details on how to reconstruct perturbation amplitudes from gauge-invariant quantities 
see [112 . 



IV. WAVE PROPAGATION ON THE ROBERTS BACKGROUND 

So, we wish to study scalar field wave propagation on the Roberts background. Typically, 
one would specify initial data for the wavepacket either on some spacelike Cauchy surface or 
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on initial null surface, and trace the later evolution using the field equations. Our choice for 
initial surface is u = const (s = 0), which forms a complete null surface if extended to the 
center of the fiat spacetime part, as illustrated in Fig. 0. The part of a pulse propagating 
through fiat background evolves trivially, and can be equivalently replaced by specifying 
field values on f = hypersurface. Thus, in our problem, the initial conditions for the 
linearized Einstein-scalar field equations are given on the s = surface, while the boundary 
conditions are determined by junction conditions across the null shell v = 0, as outlined in 
Appendix 0, and the requirement that the perturbations be bounded at future infinity. 

As was shown earlier, the Einstein-scalar field equations for the spherically symmetric 
perturbations of the Roberts solution can be reduced to the master differential equation 

Vf{x,s) = (25) 

for the single gauge-invariant quantity f{x,s) describing perturbation of the scalar field. 
The explicit form of differential operator V in scaling coordinates, given by equation (pO|), 
is 

Vf = {l- e~'n ^ + 2^ + 2^-4/. (26) 

OX^ ox OS OS 

Because of the scale-invariance of the background, the coefficients of the differential operator 
V do not depend on scale, and so the problem can be reduced to one dimension by applying a 
formal Laplace transform with respect to the scale variable s to all quantities and operators. 
In particular, for Laplace transform of / we have 

oo 

F{x,k)= J f{x,s)e-''ds, (27) 
with the inverse transformation being 

/(a:,s) = ^ / (28) 

K—ioo 

The Laplace transform can be done provided that / can be bounded by an exponential 
function of s (that is, there exist constants M, kq such that \f{x,s)\ < Me'^°^), which is 
a physically reasonable condition. The contour of integration in the complex fc-plane for 
the inverse transform ( pBD must be taken somewhere to the right of (k > > 0). The 
properties of functions of complex variables will guarantee that the result of integration is 
independent of the particular contour choice. 

When applying the Laplace transform to the differential operator. 







ds 



kF-f{s = 0), (29) 



so the initial conditions of the original problem will enter as source terms on the right hand 
side. Therefore the Laplace transform of the equation (pSf) is 
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VkF{x,k) = h{x), 



(30) 



where = CgT) is now an ordinary differential operator, algebraic in /c, and h contains 
information about the initial shape of the wavepacket at s = 0. Boundary condition on 
equation (|30| ) are inherited from the original problem by Laplace transformation. 

The explicit forms of the operator Vk and the relationship of / to the perturbation 
amplitudes are the simplest when expressed in slightly different spatial coordinate, related 
to the old one by 



l/ = e^^ = l-^. (31) 
The differential operator T>k is hypergeometric in nature 

V, = yil - y)^ + [l-ik + l)y]^ - {k/2 - 1), (32) 



with coefficients 



c = 1, a + b = k, ab = k/2 — l, 



a,b= l/2{kTVk^ -2k + 4). (33) 
The right hand side h depends on the initial conditions as 

hiy) = -yf{y,s = 0)-'^f{y,s = 0). (34) 

Here and later dot denotes derivative with respect to y (' = d /dy). Once the solution for 
F (and hence its inverse Laplace transform /) is found, one can reconstruct the other two 
gauge invariant quantities a and p describing the metric perturbations using equations (pID, 
and, in principle, write expressions for perturbations in any desired gauge choice. 

Thus, the study of wave propagation on the Roberts background is reduced to solving 



the inhomogeneous hypergeometric equation (pOf). The following analysis relies heavily on 



certain properties of the hypergeometric equation, which are collected in Appendix |^ for 
convenience. 

So far, we have not talked about specifics of the boundary conditions placed on the 
equation (^). That depends on the physical problem being considered. If the flat spacetime 
part f < were unperturbed, the null shell junction conditions, discussed in Appendix 0, 
would require that / = on the surface f = 0. If some part of the pulse propagates in the 
flat sector v < 0, f should be continuous across the surface v = 0. Essentially, we can specify 
value of / on the wedge (f = 0) U (s = 0) arbitrarily, keeping in mind that perturbation 
value should be bounded at future infinity. It is practical to split the wavepacket to three 
components, as shown in Fig. ^ and consider outgoing, "constant" and incoming packets 
separately. 

A. Outgoing Wavepacket 

The outgoing wavepacket is characterized by 
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f{x = 0,s) = fois), f{x,s = 0) = (35) 

and is propagating outwards to future infinity, except for backscatter on tlie background 
curvature which goes towards the singularity at s = +00. The boundary conditions and the 
initial term for equation (|30D are 

F{y = l,k) = Fo{k), h{y) = 0. (36) 

The general solution of the homogeneous form of equation (|30D is 

F{y, k) = A{k)Z,{y, k) + B{k)Z2{y, k), (37) 

where Zi and Z2 are two linearly independent solutions of the homogeneous hypergeometric 
equation, in notation of Appendix |B|. To satisfy boundary conditions at y = 1, parameters 
A and B must be 

A{k) = Fo{k), B{k) = 0. (38) 
Therefore, the outgoing wavepacket solution is given by 

/(z/,s) = — / Fo{k)na,b;k;l-y)e'''dk, (39) 

K—ioo 

where is the hypergeometric function. 

If fo{s) does not grow exponentially as s — > +00 by itself, i.e. the image Fo{k) does 
not have poles in Re/c > half-plane, than neither does /(y, s). The outgoing wavepacket 
just propagates harmlessly out to future infinity, never growing enough to cause significant 
deviation of the solution from the Roberts background. 

B. "Constant" Wavepacket 

Even more trivial is the case of the "constant" wavepacket, characterized by 

/(a; = 0,s) = C = /(x,s = 0). (40) 
The boundary conditions and the initial term for equation (^) are 

F{y = l,k) = C/k, h{y) = -C/2. (41) 
The general solution of equation (|30|) with these boundary conditions is 



F{y, k) = A{k)Z,{y, k) + B{k)Z2{y, k) + (42) 



and the boundary conditions at y = 1 require that 



k{k-2) 

So the constant wavepacket solution is given by 



A{k) = . , B{k) = 0. (43) 



K+JOO 



Im J 

K—ioc 

and it does not grow as s ^ +cxd either. 



2 

1 - -j^^i^^ '^-y) 



k-2 



dk, (44) 
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C. Incoming Wavepacket 

By far, the most physically interesting case is the incoming wavepacket characterized by 

/(x = 0,s) = 0, f{x,s = 0) = fi{x). (45) 

It propagates directly towards the singularity and is responsible for near-critical behavior 
and breaking of the solution away from the Roberts background, as we shall demonstrate. 
The boundary conditions and the initial term for equation (^) are 

F(x = 0,k) = 0, hiy) = -yfi{y) - fj{y)/2. (46) 

To solve the inhomogeneous hypergeometric equation (POD, 

1/(1 - y)F +[l-{k + l)y]F - {k/2 - l)F = h, (47) 

with the boundary conditions 

F{y = l,k) = 0, F{y = oo, k) bounded, (48) 

we must construct a Green's function out of the fundamental system of solutions of the 
homogeneous equation 

Ziiy) = J^{a,b; k;l-y) 

Z^{y) = (1 - yf-^'J^il - a, 1 - 6; 2 - A;; 1 - y), (49) 

where parameters a and h of hypergeometric equation depend on k as given by equation 
(^). The Wronskian of the above system is 

W{y) = {k-l)y-\l-y)-\ (50) 

and the Green's function is constructed as 

G(y, r/) = AZM + BZM ± . , . [ZMZM " zMZM] 

2po{v)W{r]} 

= AZM + BZ,{y) ± ^^^^^ [ZMZM - Z2{y)ZM] , (51) 

where the coefficients A and B are to be determined by applying the boundary conditions, 
and the plus-minus sign is taken depending on the arguments of the Green's function 

± = |+' ^^^^^ . (52) 

The Green's function G satisfies VkG{y,ri) = 6{y — r]), and hence can be used to construct 
the solution of inhomogeneous equation 

oo 

F{y,k)= J G{y,v)h{v)d7^ (53) 
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As the initial value problem (^T]) is not self-adjoint, the Green's function ( pTD need not 
be symmetric in its arguments y and rj. Note that the Green's function is calculated for 
particular fc-mode, and so depends on /c, but we omitted the third argument in 77; k) 
for brevity. 



We now proceed to apply boundary conditions to the Green's function (0), starting 
with the boundary conditions at y = 1. The fundamental solution Zi goes to one there, 
while the behavior of Z2 is fundamentally different depending on the sign of Re(l — k). If 
Refc < 1, the real part of power of (1 — y) in ( ^9]) is positive, and Z2 goes to zero when 
y = 1. If Re /c > 1, the real part of power of {1 — y) is negative, and hence Z2 diverges when 
y = 1. Substituting this into Green's function, we get 



G{y = l,r,)=A+ \^ _ ^ Z^iv) + 



0, Refc<l 
00, Rek > 1 



(54) 



The boundary conditions (|i8|) require that G{y = 1,1]) = 0, which uniquely fixes the coeffi- 
cients 

provided Refc > 1, which is precisely the region of the complex fc-plane the contour of the 
inverse Laplace transformation should be in. (If Kek < 1, coefficient B can be arbitrary.) 
With these coefficients, the Green's function (^) becomes 



[Z^{7])Z2{y)-Z2{v)Z^{y)], v<y 



The causality of wave propagation is apparent here: the wave at y is only influenced by the 
data from the past 1] < y. 

We see that the boundary conditions at y = 1 already fix Green's function (^6|), but 
we still have to satisfy the boundary conditions at infinity! One can show that they are 
satisfied automatically if (and only if) Re a > 0. Curves Refc = 1 and Rea = split 
the complex k plane into several regions, as shown in Fig. |^. The Green's function ( |56D 
as written above is defined in region A, but could be analytically continued to the whole 
complex space. The obstructions Green's function encounters on the boundaries between 
A and B and A and F are not poles, indeed, they are not even singular for regular {y,ri). 
They are rather caused by the fact that the Green's function ( pU] ) fails to be applicable 
once you cross these boundaries; in the region B boundary conditions at infinity fail to be 
satisfied, and in the region F free modes (solutions of homogeneous equation, that is) exist 
that satisfy all boundary conditions, making the Green's function not unique. The existence 
of free modes in region F is at the heart of the matter, as they grow and will determine what 
will happen to the wavepacket at the later times. (It is also possible to construct Green's 
function in region C, but since it has no bearing on our analysis, we will not do it here.) 

Once Green's function has been determined, it is simple to construct later-time evolution 
of the wavepacket from the initial data using equation (|53D 



y / -| \k—i 

^) = / k-1 - ^2(r/)^i(y)] Hv) dv. (57) 
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To get back from the complex fc-plane dependence to the physical time evolution, one per- 
forms inverse Laplace transformation 

K+ioo 

/(2/,5) = ^ / (58) 

We emphasize again that a particular choice of the contour of integration is not important, 
as long as it is to the right of the obstructions on the complex plane, in our case region F. 
In practice, one chooses the contour so that the integral (|58D is easier to evaluate. For some 
approximation to work, the contour should touch the obstruction, which means pushing it 
leftwards to the very edge of region F at Re A; = 1. 

V. LATE-TIME BEHAVIOR OF INCOMING WAVEPACKET 

While expressions for / written down in the previous section formally solve the problem 
of wave propagation on the Roberts background, they are too complicated to be of practical 
use. In this section, we use the method of the stationary phase to obtain late-time (large s) 
asymptotic for f{y, s) and analyze several physically important regimes of the wavepacket 
evolution. 

The method of the stationary phase deals with the approximate evaluation of Fourier- 
type integrals 

/(A) = J F{k) exp[iA5(A;)] dk (59) 

a 

for large positive parameter A. It is based on a simple idea that where exp[iAS'(/c)] is 
oscillating extremely rapidly and F{k) is smooth, the oscillations will cancel out, and the 
only contributions to the integral will be from stationary points of phase S{k), singular 
points of F{k) and S{k), and possibly end points. 

Inverse Laplace transform integrals (^B|) are precisely of the above type, with phase 
S{k) = k and large parameter A being the scale coordinate s. So the stationary phase 
method tells us that the asymptotic s ^ cxd of the solution f{y, s) is given by singular 
points of F{y, k) as a function of k, i.e. singular points of Gijj, rj; k). Therefore the study of 
analytic properties of Green's function (|56D plays a key role in understanding the late-time 
evolution of the wavepacket. The possible sources of non-analyticity in Green's function are 
listed below: 

• Branch points k = l± oia,b= l/2{k ± Vk^ - 2k + 4) 

• Poles at k = 2 + n in Z2 and k = —n in Zi coming from the hypergeometric function 

• The pole at A; = 1 from the prefactor 

• Power-law singularity of the type (1 — y)''^^ 
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Problem with branches of the coefficients a, b is absent in the Green's function G because 
they only enter it through the first and second arguments of the hypergeometric function, 
and the hypergeometric series are written in terms of ab = k/2 — 1 and a + b = k only. 
Various poles at integer values of k are all canceled out because of the antisymmetric way 
hypergeometric functions enter G. In fact, the only source of non-analyticity in G is power- 
law singularity, and then only ai y,r] ^ 1. Despite the appearance, G{y,ri] k) is an entire 
analytic function of k provided that y, i] are regular points. 

Since the small 1] region is important for the late-time evolution of the wavepacket, it 
is instructive to take a closer look at the approximation to the Green's function ( pBD there. 
Using the asymptotic behavior of Zi, near = 1, given in Appendix O, we obtain 



G{y,V 



k-i 



k-1 



k-1 



1 — rj 



1-y 



Z2{y)-{l-vy-'Z^iy) 

k-1 

^2(y)--Fi(y) 



(60) 



where we introduce the short-hand notation 



J^i{y) = T{aM k] 1 - y), 
J^2{y)=n^-a,l-b;2-k;l-y). 



(61) 



Note that even though poles in the hypergeometric functions no longer cancel in (|60|), they 
are purely artifacts of the approximation, and should be ignored. 

Now we will use the above approximation (|60| ) for the Green's function to study the 
late-time evolution of the packet in several important regimes. 



A. Evolution near v = 



Let us first consider the behavior of the wavepacket near y = 1, that is, near the initial 
null surface v = 0. We take the asymptotic behavior of the initial term h to be the fairly 
generic power law 



h{r])cx{l-r]r, 



(62) 



which covers the usual case of functions analytic at y = 1 (via Taylor expansion), as well as 
the case of functions with a power singularity at y = 1, such as free modes Z2. Then the 
approximation to the Green's function (|60|) gives the solution in the desired region. 



F{y,k) 



1-k 



V 



k-1 



-^2(2/)--Fi(y) 



1 — 7] 



[l-yTdT]. 



(63) 



The above integral can be explicitly evaluated using the change of variable 

1 — ri\ 



C = ln 



(64) 
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which leads to the answer 





F{y,k) ^ - 



1-k 



-(1-1/) 



l+a 



k-l 



k + a l + a 



(65) 
(66) 



To obtain the late-time evolution of the wavepacket near y = 1, we use the method of 
stationary phase. Observe that the only real singularity of F{y, k) is the simple pole at 
k = —a, with the rest being artifacts of the approximation (0). Therefore, the late-time 
behavior of the wavepacket is exponential in the scale s; indeed, it is given by 



f{y,s) 



;i-i/) 



l+a 



1 + a 



■^2(1/; -a)e-"^ 



1 



l + a 



^2(2/; -a) e" 



(67) 



There are several things worth noting about the above result. First, it gives the correct 
answer for the evolution of the free mode. Since ^2(1/;^) ~ (1 — y)^"^, the initial term works 
out to be h{y) ~ (1 — p)(l — y)~^, so the above formula gives f{y, s) ~ ^2(1/; p)e^'^, which is 
precisely what the evolution of the free mode should be. Second, the growing modes cannot 
be excited by the initial profile analytic at ?/ = 1. If this were the case, then h could be 
expanded in a Taylor series around y = 1, each term in the series raised to integer power 
corresponding to non-negative value of a, and so not growing at large s. Third, the above 
argument illustrates how only the Z2 content of the initial data is relevant to the subsequent 
evolution of the wavepacket. 

So it seems that the exponential growth of the wavepacket at late times must be already 
built into the initial data in form of a power-law divergence of the initial wave profile, just 
as it is encoded in the pure free mode Z2. But a power-law divergence of the perturbation 
near v = causes curvature invariants to diverge at the junction, making the surface v = 
a weak null singularity, which casts a shadow of doubt on the physicality of such growing 
modes. The question then arises whether it is possible to somehow eliminate the offending 
divergence, while still having a wavepacket grow to large enough values. The answer to this 
question is yes, and it is discussed next. 



B. Evolution of a wavepacket initially localized at v = 

So what will happen if we cut off a diverging initial wave shape (|62|) below some small 
value of y — 1, say A? To put it in another way, will the power-law diverging wave localized 
near y = 1 backscatter and affect the evolution of the wavepacket at the large y? If not, 
then subtracting such localized wavepacket from perturbation modes discussed above will 
cut off the divergence at y = 1, while keeping the rest of the wavepacket evolution essentially 
unchanged. 

We model such a localized wave by adding an exponential cutoff to the generic power 
law initial term ( p^ ) 

h{r]) oc {l-r])"e'^^-'^^/\ (68) 
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The exponential factor is chosen because it effectively suppresses h for values of y—1 > A, yet 
still keeps the calculations simple. We are interested in the evolution of the wavepacket well 
outside the region of initial localization, but still for small enough y so that approximation 
(|60|) holds, that is for A <^ ?/ — 1 <^ 1 (that can always be arranged for small enough A). In 
this region of interest we have 



F{y,k) 



k 



1—1]" 



k-l 



T2{y)-T,{y) 



a(l-ri)/\ 



drj. 



The above integral can be evaluated by the change of variable 

1] — 1 



t 



A 



which yields the following approximation for F{y; k) in the region of interest: 

I ^2(2/)--^l(y) 



F{y,k) 



1 



1 - k 
1 



-At 



At)"e-* Xdt 



Xf^^Tik + a){l- yY-':FM - (-A)^+"r(l + a)^i(y) 



(69) 



(70) 



(71) 



(72) 



So outside the region of initial localization of the wavepacket, but still for small the 
Laplace transform of the field perturbation is approximately given by 



F{y,k) 



1 



k-l 



\f^^T{k + a)ZM - (-A)^+"r(l + a)Z,{y) 



(73) 



The main contribution to the late-time behavior is coming from the poles of the gamma- 
function T{k + a) in the first term. Using the stationary phase approximation, it is possible 
to calculate this contribution exactly. The inverse Laplace transform of F{y\ k) can be 
reduced to the inverse Mellin transform of the gamma-function 



Z2iy; -a) e" 



a 

ft+a+ioo 



1 + a 27ri 
Z2{y;-a) e— 



K+a~ioo 
ft+a+ioo 



1 + a 2ni 



\^V{k)e^'dk 

m 



(s+lnA) 



-k 



dk 



K+a—ioo 



Z2iy; -a) 

l + a 



-(s+lnA) 



(74) 

(75) 

(76) 
(77) 



which is a known integral; indeed, [M. ^r](a;) = e Therefore we obtain the following 
late-time approximation of the field perturbation outside the region of initial localization: 



e(s + lnA) 

l + a 



^2(2/; -«) e" 



(78) 
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This looks very similar to the earlier result ([67|); the only difference is the factor 0(s + lnA), 
where is defined by 

e(x) = exp [-e-^] , (79) 

and is almost a step function, rapidly changing its value from to 1 as its argument becomes 
positive 

where the width of the transition is of order unity. This means that the perturbation outside 
the region of initial localization does not feel the effect of the field at ?/ — 1 < A until a much 
later time, namely s = — InA, when it suddenly spreads. To put it simply, the wavepacket 
initially localized at v < X does not backscatter until it hits the singularity at m = 0, and 
then goes out in a narrow band —u < A. 

The above result shows how we can cut off the free mode Z2{y; k) to avoid the curvature 
divergence, and yet have it grow sufficiently large. To quantify how large can it grow, consider 
the free mode / = Z2{y] fc)e'^* (1 — yY~^e^^, with divergent scalar curvature R oc /', and 
cut it off at ?/ — 1 < A. The largest initial curvature value is of order R oc A~'^, while the 
initial energy of the pulse is \R oc A^"*^, and can be made arbitrarily small. The perturbation 
mode will grow exponentially until the cutoff backscatters at s = — In A, at which time its 
amplitude will be A~*^, with proportionally large curvatures and energies. In other words, 
the initial large curvature seed localized in a small region spreads over the whole space in 
the course of the evolution, with the energy of the pulse growing correspondingly. Thus, the 
free perturbation modes considered here are physical and grow exponentially to very large 
amplitudes, certainly enough to leave the linear regime, and are therefore responsible for 
the evolution of the solution away from the Roberts one. 



C. Generic initial conditions 

We now turn our attention to the evolution of the wavepacket from generic initial condi- 
tions. It is reasonable to expect that completely generic initial conditions will have non-zero 
content of all perturbation modes present in the system, both growing and decaying, as 
given by equations ( p7D and 

f{y.s)= J [W2{k)Z2{y;k) + W,{k)Z,{y;k)]e'^dk. (81) 

K—ioo 

However, decaying modes will disappear very quickly, so only the growing modes are relevant 
to late-time evolution. Assuming that the content of the free growing mode ^2(1/; k) is given 
by the weight function W{k), the generic wavepacket evolution is given by the sum of all 
such modes 

fiy,s)= J Wik)Z2iy;k)e''dk, (82) 
r 
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where the infinite contour of integration F runs vertically in the regions F and F of the 
complex plane on Fig. at Re/c = 1. However, we note that the part of the contour 
between endpoints = of the regions F and F does not correspond to free growing 

modes, as Z2{y — > oo) fa (— with Rea < there, and so boundary conditions at infinity 
are not satisfied. Therefore, for the initial wavepacket bounded at infinity, the content of 
such modes is suppressed, so that / dk [W{k) {—y)~'^] ~ 1 for large y. This leads to slower 
growth rates / dk [W{k) (—|/e*)~"e^'^ ~ e^^^"-** at the later times. Hence, the piece of the 
contour F between the endpoints k^ can be omitted from the integration without affecting 
the late-time evolution. 

For completely generic initial conditions, we should expect to be a smooth function 
of k in the free mode region F, not preferring any particular value of k. Therefore, using 
the stationary phase approximation, the main contribution to the late-time behavior of the 
above integral comes from the end points of the contour of integration, 

/(y, s) ^ -W{k+)Z,{y- k^) ^ - W{k^)Z,{y- k^) (83) 

Ignoring the overall weight factor, we find that the late-time evolution of the generic 
wavepacket is given by 



f{y,s) oc Re 



Z2iy;ko] 



gfcos 



(84) 



We emphasize that the single fco-mode dominates the course of evolution of the generic 
wavepacket, and thus a certain universality is present in the way a generic perturbation 
departs from the Roberts solution. 



VI. DISCUSSION 

In this paper, we studied spherically symmetric perturbations of the Roberts solution 
with the intent to understand how nearby solutions depart from the Roberts one in the course 
of the field evolution, and what bearing the Roberts solution has on the subject of critical 
phenomena in general, and how it is related to Choptuik's critical solution in particular. We 
analyzed the behavior of incoming and outgoing wavepackets, and we focused our attention 
on the incoming one as the physically relevant one for the question posed. With the aid of 
the Green's function formulation, we were able to completely solve the perturbation problem 
in closed form, as well as obtain simple approximations for the late-time evolution of the 
field in several important regimes. 

As was shown above, the departure of the generic perturbation away from the Roberts 
solution is universal in a sense that the single mode Z2{y] fco)e^°* dominates the late-time 
evolution of the field. The complex growth exponent gives rise to an interesting physi- 
cal effect: the perturbation developing on the scale-invariant background evolves to have a 
scale-dependent structure cos(lm kos) . The exponential growth of the amplitude of the 
perturbation will eventually be stopped by the non-linear effects, while the periodic depen- 
dence of the perturbation on the scale will most likely remain. The period of oscillation, 
obtained in the linear approximation, is 
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What does this periodic dependence of the solution on the scale mean physically? To 
answer this question, lets see how this symmetry is expressed in the Schwarzschild coordi- 
nates {r,t) often used in numerical calculations (see, for example, In this coordinates 
the metric (||) is diagonal 

ds"^ = -adt^ + pdr"^ + r^dn\ (86) 

where metric coefficients and explicit expressions for coordinates are given in Appendix 0. 
For our purposes, it suffices to note that the coordinate x determines the ratio r/t, while 
the coordinate s sets an overall scale of both space and time coordinates via e"'* factor 

r = exp[x — s]. (87) 

One can see that taking a step A in the scale variable s is equivalent to scaling both spatial 
and time coordinates r and t down by a factor e^. Therefore, the solution being periodic 
in scale coordinate s is equivalent to being invariant under rescaling of space and time 
coordinates r and t by a certain factor 

fix, s + A) = fix, s) ^ fie-\ e-H) = fir, t). (88) 

The later is an expression of the symmetry observed in the numerical simulations of the 
massless scalar field collapse |^ , and referred to as echoing, or discrete self-similarity in the 
hterature 

Thus, our simple analytical model of the critical collapse of the massless scalar field 
illustrates how the continuous self-similarity of the Roberts solution is dynamically broken 
to discrete self-similarity by the growing perturbations, reproducing the essential feature of 
numerical critical solutions. The value of exponent for the endpoints of the spectrum of 
growing perturbation modes, kQ = 1 + i\/2, gives the period of discrete self-similarity as 

A = V27T = 4.44 (89) 

for linear perturbations of the Roberts solution, which is within 25% of the numerical value 
A = 3.44 measured by Choptuik Given that the perturbative model considered here 
reproduces all the symmetries of the Choptuik's solution, and gives a good estimate for 
the period of echoing, it is instructive to compare the actual field profiles to the numerical 
calculations. There are some technical issues connected with rewriting our results in the 
variables Choptuik uses, which are addressed in Appendix 0, but the end result of calculation 
for field variable X = \/2tt |^ from the perturbation modes is presented in Fig. |^. 
Comparing this plot to Fig. 2 in original Choptuik's paper we see that they share 
one common feature, which is the oscillatory nature of the field solution; however, the 
shape of the field profiles is quite different. This discrepancy is not surprising, however, 
since perturbation methods in critical phenomena are usually viable for calculating critical 
exponents, but not the field configurations themselves. 



r 

1 



exp 



x + -e 



2x 
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The emerging discretely self-similar structure, and the universal way in which the generic 
perturbation departs away from the Roberts solution offer support for the conjecture that 
the Roberts solution is "close" to the Choptuik one in the phase space of all massless scalar 
field configurations (in a sense of being in the basin of attraction of the latter), and will 
evolve towards it when perturbed |T^. It seems highly unlikely, however, that the critical 



mode responsible for the decay of the Choptuik solution will be completely absent in the 
initial data originating near the Roberts solution, as this usually requires fine-tuning of the 
parameters. So though at first the field configuration near the Roberts solution might seem 
to evolve towards the Choptuik solution, after a while the critical mode will kick in and drive 
the field to either dispersal or black hole formation. This picture is in line with the Choptuik 
solution being intermediate attractor, and we offered "ball rolling down the stairs" analogy 
earlier |]T^ to visualize the field evolution as it goes from initial configuration (the Roberts 
solution) to local attractor (the Choptuik solution) and then to global attractor (black hole 
or fiat spacetime) in the phase space. Unfortunately, linear perturbation methods are not 
sufficient to provide a proof of the proposed scenario, and fail to give the answer as to what 
would the eventual fate of the evolution be (whether black hole or flat spacetime end-state 
will be selected), and how fast would the field get there. 

To completely answer these questions, one would need to employ some sort of non-linear 
calculation, or perform numerical simulations of the evolution. In particular, it would be 
interesting to evolve the perturbed Roberts spacetime numerically and look for the Choptuik 
spacetime as the possible intermediate attractor. Nevertheless, it may happen that some 
information about Choptuik's solution can be gained from the linear perturbation analysis 
of the Roberts solution. The appeal of this method lies in the fact that such analysis could 
be carried out analytically, while Choptuik's solution is still unknown in the closed form. 
Similarly, one can try to study properties of other analytically-unknown critical solutions in 
different matter models based on "nearby" solutions with higher symmetry and simpler form. 
One might also hope to obtain acceptable analytical approximations to critical solutions, 
Choptuik's in particular, by going to higher order perturbation theory in the region near 
the singularity. 

To summarize, our main result is the dynamical explanation of how discretely self-similar 
structure forms on the continuously self-similar background in the collapse of the minimally- 
coupled massless scalar field, and the prediction for the period of this structure A = y^ir = 
AAA, which is quite close to the numerical value A = 3.44, considering we only did a first- 
order perturbation analysis, and that non-linear effects can, in principle, renormalize that 
value. 
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APPENDIX A: NULL SHELL JUNCTION CONDITIONS 



The influx of the scalar field in the Roberts solution is turned on at the advanced time 
V = 0, with the spacetime being flat before that. In this appendix, we consider in detail 
junction condition across this null surface, and their implications for boundary conditions of 
perturbation problem. This discussion is adopted from general treatment of thin null shells 
by Barrabes and Israel . 



Consider the general spherically symmetric metric in null coordinates 

ds^ = -26^" dudv + r^dn^. (Al) 

To fix the geometry of the soldering of spherically symmetric spacetimes uniquely, one must 
match radial functions across the constant advanced time hypersurface. To this end, we 
rewrite the above metric in terms of advanced Eddington coordinates on both sides 

ds^ = -e^dv{fe^dv-2dr) + r^dQ'^. (A2) 

The metric coefficients in Eddington coordinates can be calculated from those in null coor- 
dinates by 

e^ = - — , / = -2^. (A3) 

The surface density and pressure of the null shell are then determined by jumps of the metric 
coefficients across the v = const surface 

inr^e = [m], SnP = [ip^r], (A4) 

where the local mass function m{v,r) is introduced, as usual, by 

/=l-^. (A5) 

r 

For Minkowski spacetime / = 1 and ip = const, so below the v = hypersurface surface we 
have 

m{v < 0) = 0, ip^riv < 0) = 0. (A6) 

For the Roberts solution we have a = 0, = u"^ — uv, so a direct calculation gives 

uv . 1 v"^ , , 

= ~iry 'V,r = - -r^, — (A7) 



Obviously, 



lim m = 0, lim -0^ = 0, (A8) 



SO the Roberts solution is indeed attached smoothly to the flat spacetime, without a delta 
function-like stress-energy tensor singularity associated with a massive null shell. 
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If one wishes to attach a perturbed Roberts spacetime to the flat one, as we do, and 
still have no singularity at the junction, the null shell matching conditions above will place 
boundary conditions on the perturbation values at the junction surface. To obtain these, 
we take the perturbed Roberts solution in null gauge of Ref. |T^, given by a = a{x)e''^, 
r = e^^^{l + p{x)e^'^), and calculate the null shell surface density 

[m] = [-p'/2 + ik/2 - l)p + crJeC^-^)^ (A9) 

and surface pressure 

[^^,] = -k\{k - l)p + 2(T]e('=+^)^ (AlO) 

The simultaneous vanishing of these two for arbitrary k can only be accomplished if 

p = p' = = 0, (All) 

and these are precisely the boundary conditions we imposed on metric coefficients earlier. 
One also has to require continuity of the scalar field across junction, so the boundary con- 
dition on scalar field perturbation is 

= 0. (A12) 

These boundary conditions on perturbations in the null gauge simply mean vanishing of 
gauge-invariant perturbation amplitudes (defined in the next appendix) on the f = hyper- 
surface. 



APPENDIX B: PROPERTIES OF THE HYPERGEOMETRIC EQUATION 

We showed above that the linear perturbation analysis of the Roberts solution can be 
reduced to the study of solutions of the hypergeometric equation with certain parameters. 
The hypergeometric equation has been extensively studied; for a complete description of 
its main properties see, for example, [|T6|. In this appendix we collect the facts about the 
hypergeometric equations that are of immediate use to us, mainly to establish notation. 

The hypergeometric equation is a second order linear ordinary differential equation, 

y{\ - ?/)Z + [c - (a + 6 + \)y\Z - ahZ = 0, (Bl) 

with parameters a, b, and c being arbitrary complex numbers. It has three singular points 
at y = 0,1,00. Its general solution is a linear combination of any two different solutions 
from the set 

Zi = jF(a, b; a + b + 1 — c; 1 — y) , 
^2 = (1 - yy-'^'^'J^ic -a,c-b;c+l-a-b;l-y), 
Z-i = (-l/)-"^(a, a + l-c-a + l-b- y-'), 
^4 = i-y)-'Hb + 1 - c, 6; 6 + 1 - a; y'^), 
= J^(a,6; c;y), 

Ze = y'-'^Fia + l-c,b+l-c-2-c-y), (B2) 
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where T{a^ b; c; y) is the hypergeometric function, defined by the power series 



J^(a, b; c;y) = 

n=0 



(B3) 



and we used shorthand notation (a)„ = r(a + n)/r(a). The hypergeometric series is regular 
at y = 0, its value there is JF(a, 6; c; 0) = 1, and it is absolutely convergent for \y\ < 1. 
Considering the hypergeometric series as a function of its parameters, one can show that 
J-'{a,b;c;yo)/T{c) is entire analytical function of a, b, and c, provided that \yo\ < 1. 

The solutions Zi, . . . , Zq are based around different singular points of the hypergeometric 
equation, with asymptotics given by 



Zs-- 



1, ^2 = (1 - 

[-y)-", z, = 

Z5 = 1; Zq 



yY °- ^ near ?/ = 1, 
h 



-y) near y 
y^~^ near y 



00, 
0. 



(B4) 



Any three of the functions Zi, 
In particular. 



Zq are linearly dependent with constant coefficients. 



(B5) 



' Zi' 




Cl3 


Ci4 




' z^' 


.^2. 




C23 


C24 







where the coefficient matrix is given by 

T{a+b+l-c)T(b~a) 
T(b+l-c)T{b) 
T{c+l-a-b)T{b-a) - 
T{l-a)T{c-a) ^ 



Cl3 


Ci4 




_C23 


C24 _ 





3-i7ra r(a+&+l-c)r(a-fc) jTrb 

r(a+l-c)r(a) ^ 
■iTT{c-b) r{c+l-a-b)r{a-b) iT,{^c-a) 

r(i-fe)r(c-fe) ^ 



(B6) 



These relationships are true for all values of the parameters for which the gamma-function 
terms in the numerators are finite, and all values of y for which corresponding series converge, 
with Imy > 0. If Imy < 0, signs of arguments in the exponential multipliers should be 
inverted. We shall not give the rest of similar relationships here. 



APPENDIX C: SPACETIME IN CURVATURE COORDINATES 



For comparison of our results to Choptuik's numerical simulations [jT|, we must rewrite 
them in the diagonal Schwarzschild-like coordinates Choptuik uses 

ds^ = -adt^ + fidr'^ + r'^ dn^. (CI) 

One can straightforwardly check that the coordinate change 

r = exp[2; — s] (C2) 

diagonalizes the Roberts metric (^. By self-similarity, the quantity t/r, as well as the metric 
coefficients a and (3 do not depend on the scale s, but only on the coordinate x. The metric 
coefficients, written as functions of x, are 



t = — exp 



-s e 

2 



2x 
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a = 2^^, /3 = 2 \ . (C3) 



If one wishes, one can rewrite them as exphcit functions of t/r, using 

X = ^\nW{ryt^), (C4) 

in terms of Lambert's VT- function, which is defined by the solution of transcendental equation 

Wexp{W)=x. (C5) 
The expressions for metric coefficients are then 

exp[WV^ . o Wjryt') 

However, the coefficients a and f3 cannot be written in closed form in terms of elementary 
functions of t/r. 

As you can see from expressions for the metric above, diagonal Schwarzschild coordinates 
are not particularly well-suited for description of the Roberts spacetime. On top of the 
complicated metric form, one artifact of the diagonal coordinate system is that the null 
singularity at m = gets compressed into a point at r = t = 0. Also, slices t = const cut 
across the f = hypersurface, so one has to be careful with discontinuities of the solution 
there. 

The perturbation amplitudes in the gauge preserving diagonal form of the metric are 
also quite complicated. The simplest way to get them from gauge-invariant quantities is to 
explicitly a find gauge transformation 

e^ = (Ai?,0,0) (C7) 

connecting the simple field gauge K = k^^ = with the diagonal gauge, fixed by conditions 
K = and (2m— f )^/c^^ = u'^k^u- The effects of the gauge transformation on the perturbation 
amplitudes were given above in Section pTj Imposing the condition K = 0, one finds that 
B must be related to A by 

5 = ^A (C8) 
u 

A is then found by imposing the other condition fixing diagonal gauge, which leads to the 
following equation 

(2m - vfA^^ - u{2u - v)A^u -vA = 2uJf dv. (C9) 
Rewriting A in scaling coordinates, 

A{y,s)=A{y)e'^'-'^', (CIO) 
transforms the above equation into the ordinary differential equation 
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(l + y)i+[l-A;/2(l + y)]^=- I Fdy, 



(Cll) 



which can be easily solved to give 



A{y) 



l + y 



(C12) 



1 1 



Once the connecting gauge transformation is known, it is trivial to obtain the perturbation 
amplitudes in the Schwarzschild diagonal gauge. In particular, the scalar field perturbation 
is given by 



We end this section by observing that while the gauge transformation term A is a small 
correction to gauge- invariant quantities near y = 1, it is not at all well behaved at infinity. 
Indeed, it blows up exponentially as ea^/l The presence of this gauge artifact in the quite 
sensibly-looking diagonal gauge illustrates just how easily one can get into trouble if one is 
not working in a gauge-invariant formalism. 



F{y;k)+A{y;k). 



(C13) 
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FIGURES 




(null coordinates) (scaling coordinates) 



FIG. 1. Global structure of the Roberts solution: The scalar field influx is turned on at w = 0; 
spacetime is flat before that. The field evolution occurs in the shaded region of the diagram, and 
there is a null singularity in the center of the spacetime. 
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outgoing 



"constant" 



incoming 



FIG. 2. Wave propagation on the Roberts background: Initial conditions can be equivalently 
specified on the surface s = extending to the center of the flat part of spacetime (r = 0), or on 
the (x = 0) U (s = 0) wedge. By linearity, the wavepacket can be decomposed into three modes: 
outgoing, "constant", and incoming. 



Im k 




FIG. 3. Complex perturbation spectrum. Values of k to the left of the solid line are prohibited 
by the boundary conditions at infinity, to the right of the broken line by the initial conditions at 
y = 1. Values in the region of intersection (the shaded regions F and F) are allowed, and constitute 
the perturbation spectrum. 
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fa) k 



+ 3i 





(h) k=l + iV2 



FIG. 4. Field perturbation profiles on a slice u = const: (a) for a typical value of k inside 
region F, (b) for a value of k at the endpoint of region F. The horizontal coordinate on the plots 
is Inv. 



FIG. 5. Profile for the field variable X 

^ks 




on the slice t 



const for the dominant 



mode /(y, s) = Z2{y; k^) e / s. Compare this plot to Fig. 2 of Choptuik's paper Q. 
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